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Abstract 

Using data collected by the FOCUS experiment at Fermilab, we present a new 
measurement of the charm semileptonic branching ratio BR ( £> ~*_f_J* ^ , From a 

sample of 320 ± 44 and 11372 ± 161 D + -> p°n + v and D+ -> K--k + h + v events 

respectively, we find B *^ + ~ffi + + ) ) = °- 041 ± °- 006 ( stat ) ± °- 004 ( s y st )- 



1 Introduction 



Semileptonic decays provide an ideal environment for the study of hadronic 
matrix elements affecting the weak mixing angles from heavy flavor decays 
since the weak part of the current can be separated from the hadronic part 
of the current. The hadronic current, described by form factors, can be cal- 
culated by different theoretical methods, e.g., Lattice QCD (LQCD), Quark 
Model (QM), Sum Rules (SR), without the added complication of significant 
final state interactions. While many theoretical models predict the value for 
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the branching ratio BR(£>+->g°i+«/) [1,2,3,4,5,6,7,8,9,10,11,121, only a few ex- 

to BR(D+->K 1 '''''''' ' ' ' J' J 

perimental measurements of this ratio have been made. Furthermore, even 
though previous measurements of BR ( D+ Jf + u ) h ave su ff e red from a lack of 

statistics [13,14,15], most theoretical predictions still differ by at least 2a from 
the world average [16] for the muonic mode. In this paper, we present a new 
measurement of this branching ratio based on 320 ±44 D + — > p°fi + u events. 2 



2 Event Selection 



The data for this analysis were collected with the FOCUS experiment during 
the 1996-97 fixed target run at Fermilab. The FOCUS experiment utilized 
an upgraded version of the forward multi-particle spectrometer used by ex- 
periment E687 [17] to study charmed particles produced by the interaction 
of high energy photons, with an average energy of ~ 180 GeV, with a seg- 
mented BeO target. Precise vertex determination was made possible by two 
sets of silicon strips detectors consisting of two pairs of planes interleaved with 
the target segments [18] and four sets of planes downstream of the target re- 
gion, each with three views. Five sets of proportional wire chambers combined 
with two oppositely polarized analysis magnets completed the tracking and 
momentum measurement system. Charged hadron identification was provided 
with three threshold multi-cell Cerenkov counters capable of separating kaons 
from pions up to 60 GeV/c [19]. Muons were identified in the "inner" muon 
system, located at the end of the spectrometer, which consisted of six arrays 
of scintillation counters subtending approximately ±45 mrad [20]. 

The D + — > p°fi + u events are selected by requiring two oppositely charged 
pions and a muon to form a good decay vertex with confidence level exceeding 
5%. Tracks not used in the decay vertex are used to form candidate produc- 
tion vertices. Of these candidates, the vertex with the highest multiplicity is 
selected as the production vertex; ties are broken by selecting the most up- 
stream vertex. This production vertex is required to have a confidence level 
greater than 1% and be inside the target material. 

Muon tracks are required to have a minimum momentum of 10 GeV/c and 
must have hits in at least five of the six planes comprising the inner muon 
system. These hits must be consistent with the muon track hypothesis with 
a confidence level exceeding 1%. In order to reduce contamination from in- 
flight decays of pions and kaons within the spectrometer, we require the muon 
tracks to have a confidence level greater than 1% under the hypothesis that 
the trajectory is consistent through the two analysis magnets. 

The Cerenkov algorithm used for particle identification returns the negative 
2 Unless stated otherwise, charge conjugation is implied. 
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log-likelihood for a given track to be either an electron, pion, kaon, or proton. 
To identify the two pion candidates we require that no other hypothesis is 
favored over the pion hypothesis by more than five units of log-likelihood for 
each track. Furthermore, we require that the pion hypothesis for the track 
with charge opposite to the muon be favored over the kaon hypothesis by at 
least five units of log-likelihood. This very stringent is cut used to suppress 
background from the Cabibbo favored decay D + — >■ K~7T + fi + v where the 
kaon is misidentified as a pion. Additionally, the pion hypothesis must be 
favored over the kaon hypothesis by more than one unit of log-likelihood for 
the remaining track. 



To suppress short-lived backgrounds, we require the decay vertex to be sep- 
arated from the production vertex by at least 15 times the calculated un- 
certainty on the separation <jl and be outside of the target material by la. 
Because the target region has embedded detectors, the decay vertex is also re- 
quired to be outside the detector material. These cuts are especially effective 
at removing non-charm backgrounds, making the contribution from minimum- 
bias events negligible. Contamination from higher multiplicity charm events 
is reduced by isolating the rnrfi vertex from other tracks in the event (not 
including tracks in the primary). We require that the maximum confidence 
level for any other track to form a vertex with the candidate to be less than 
1%. 



Background from D* + — > D°tt + — > (7r~ fi + u)7i + , where the soft pion is erro- 
neously assigned to the decay vertex, is reduced by requiring M(it + it~ n + ) — 
M(n~fi + ) > 0.20 GeV/c 2 . Background decay modes with a neutral hadron in 
the final state, such as Df — > i]'fi + u — > (r]7T + 7r + )/i + z/, are reduced by requiring 
a visible mass cut of 1.2 GeV/c 2 < M(ir + ir~fi + ) < 1.8 GeV/c 2 . This cut elim- 
inates ~ 40% of the D + — >■ p°fi + u signal, but also proves to be very effective 
at rejecting background consisting mostly of kaons and pions misidentified as 
muons. 



In order to reduce systematic errors common to both modes, the D + — > 
K~7T + fi + i> events used for normalization had the same vertex and muon iden- 
tification cuts as in the D + — > p°fi + u analysis. The kaon is identified by 
requiring that the kaon hypothesis is favored over the pion hypothesis by at 
least two units of log-likelihood. The pion candidate is identified by requiring 
that the pion hypothesis for this track be favored over the kaon hypothe- 
sis. Background from D* + — > D°tt + — > (K~ fi + u)7i + is reduced by requir- 
ing M(K-7i + n + ) - M(K~fi + ) > 0.20 GeV/c 2 . A cut on the visible mass of 
1.0 GeV/c 2 < M(K~n + p + ) < 1.8 GeV/c 2 suppresses background from muon 
misidentification. 
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3 BR Determination 



The D + — > p°fi + u yield is estimated using a binned maximum log-likelihood 
fit of the 7r + 7r" invariant mass. The likelihood is defined as: 

#bins Si -m 

C = I] -^T~ x P^alty (1) 

i=i s *- 

where Sj is the number of events in bin i of the data histogram, n-i is the 
number of events in bin i of the fit histogram, and a penalty term, described 
below, is used to set a loose constraint on a known branching ratio. 

The fit histogram is composed of: binned, normalized shapes of signal and 
background components obtained from Monte Carlo simulations, S x , the num- 
ber of events as estimated by the fit for each shape, Y x , and the number of 
events that occur due to feed-down from the Cabibbo-favored decay D + — > 
K~n + p + v. The number of events in each bin is then: 



n-i =Y D +^ p o tl + u S p t) fl + l/ + ECY D +^x-7r+ M +^ e (-^' 7r / iZ/ — * PP^^k-k^v + 

Yd+->k°h+u^k°h+v + Y D +^ UJ ^+ v S bJ n+ v + Y D +S D + + Y C S C + Y M S M 

(2) 

where the terms in Eq. 2 are explained in detail below. 

As mentioned before, the shapes for the signal and background are obtained 
via Monte Carlo simulation. The Monte Carlo is based on Pythia 6.127 [21] and 
contains all known charm decays with their corresponding branching ratios and 
careful simulation of known secondary processes. After an event is generated, 
it is passed through a simulation of the FOCUS spectrometer. The events are 
then selected in the same way as the data. 

Y D +_, p Op+ u is the yield of the D + — > p°/x +z/ signal. ECY D +_ tK - 7T +p+ u is the 
efficiency-corrected yield (ECY) for D + — > K~tt + This quantity is the 
estimated number of D + — > K~ir + [j, + b' events produced by FOCUS. This, 
along with the Monte Carlo efficiency for a D + — > K~n + p + v event to be 
misidentified as a D + — > p°/i +z/ event, e(K-Kfiv — > ppu), provide an estimate 
of the amount of feed-down of this mode into our signal. The ECY is fixed 
in the fit to the value obtained from the D + — > K~n + p + v analysis used 
for the normalization mode. Y D +^ K Op +l/ is the yield of a small Kg — > 7r + 7r~ 
component. 

Yd+^h+v is the yield of D + — > c<j/x + z/, where the u could decay either to 
7T + 7r~7r or to 7r + 7r~. We use the recent CLEO-c collaboration measurements 
of the absolute branching ratio of D + — > K e + v and D + — > Lue + v [22] to set 
a loose constraint on the yield of D + — > uop + v. To this end, we add a penalty 
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term to the likelihood of the form 3 



cxp 



-\R /7? *oECY 
2\ W I K 



- Y D+ 

— >i 



(3) 



where R ;T 7»o 




error used in Eq. 3 is based on 



the errors in the branching fractions reported by CLEO-c with statistical and 
systematic errors added in quadrature to the error in the efficiency corrected 
yield for D + — > K p + v. 

Y D + is the combined yield of the modes Df — > tj'/jLu, — > rj/jtu, and Df — > 
4>jj, + u with i]' decaying to either p°7 or to r]n + n^ , r] decaying to either n + n^n 
or to 7r + 7r _ 7, and decaying to pre. These modes are generated simultaneously 
with their corresponding relative branching ratios, and a single shape, S D +, is 
obtained. 

Ye is the number of combinatorial background events where at least one of the 
three charged tracks forming the decay vertex does not belong to the vertex. 
After applying our selection criteria, this background is dominated by charm 
decays. In order to generate Sc , the combinatorial background shape, we use 
a large Monte Carlo sample which simulates all known charm decays, where 
after an event is selected the reconstructed tracks are matched against the 
generated tracks. If one of the reconstructed tracks does not belong to the 
generated decay vertex, the event is flagged as a combinatorial background 
event. 

The last term of Eq. 2, Y M) is the number of events due to muon misidentifi- 
cation. The muon misidentification shape is also obtained from a large Monte 
Carlo sample where all known charm decays are simulated. In this case tracks 
within the acceptance of the inner muon system with a confidence level less 
than 1% and momentum greater than 10 GeV/c are taken as muons. This 
allows us to use the same selection as in the analysis, but with very few 
real semi-muonic decays in the sample. This shape is then weighted with a 
momentum-dependent misidentification probability function to obtain the fi- 
nal shape used in the fit. The same technique, applied to a sub-sample of the 
FOCUS data, gave a shape in very good agreement with the shape used in 
the fit. We choose not to use shape or Y M estimates from the data due to the 
limited statistics available. Instead, we allow Y M to float freely in the fit. 

The D + — > K~7i + fi + u yield used for the normalization is estimated in the same 
way as the D + — > p°fi + u yield. In this case, we only need two components in 
the fit: one for the signal and one for the background. This background shape 
is obtained from a Monte Carlo sample where all known charm decays, except 
D + — > K~ir + p + is, are generated. The fit to the K~ir + line shape is similar 

3 Here we have assumed that the electronic and muonic rates are equal. 
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to the one used in [23], and though this fit does not describe the complex 
line shape presented in [24], we find that it provides a robust estimate of the 
D + — >• K~it + n + v yield. 

From our fit we find 320 ± 44 D + -> p°fi + u events and 11, 372 ± 161 D+ -> 
K~n + fi, + is events. In Fig. 1 we show the fit result. The ratio of branching 
ratios is defined as 

BR(£>+ -» pV+z/) = Yp+^^/ep+^ffi^+u ^ B R(X* -> 7^+) 

— »K TT + fl+U 

(4) 

This branching ratio must be corrected to account for the (5.30 ± 0.74lg;gg)% 
non-resonant S-wave contribution present in the D + — > K~7r + fi + v spectrum 
previously reported by FOCUS [20,24]. Including this correction, with errors 
added in quadrature, we find 

BR(D+ -> pV + ") _ . 0412±0 . 0057 . 



BR(D+ -> K"°ii+v) 



Decay Mode 


Total Yield 


Yield in signal region 


D + -► pV + ^ 


320 ± 44 


282 




68 a 


44 


D+ -► irg/i+i/ 


7±6 





D+ modes total 


179 ± 40 


101 


D+ -► 


51 ±22 


10 


Muon Mis-Id 


550 ± 44 


263 


Combinatoric 


233 ± 50 


99 



a The D + — > K~n + yield is not a fit parameter, instead it is estimated based 
on the efficiency for a D + —> K~ir + n + v event to be reconstructed as a D + —> 
p°p + u event as described in the text. 
Table 1 

Contributions to the tt + tt~ invariant mass spectrum. The third col- 
umn shows the number of events in the signal region defined as 
0.62 GeV/c 2 < M{-k+h~) < 0.92 GeV/c 2 . 

Many tests were performed to ensure that the final result is stable to our 
detailed cut choice, as well as a good representation of the data. The quanti- 
tative tests we performed to determine additional sources of uncertainty are 
detailed below. 
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M(lt + 7I ) GeV/c ! M(K"lt + ) GeV/c ! 

Fig. 1. Fit to the 7r + 7r~ invariant mass (left) and the K~tt + invariant mass (right). 
The fit to the data (error bars) is shown as a solid line. The hatched histogram on 
invariant mass plots is the background. 




M(lt + 7I ) GeV/c ! M(lt + 7I ) GeV/c ! 

Fig. 2. M(7r + 7r~) background contributions shown in cumulative plots. Left: semilep- 
tonic contributions. Right: Muon misid and combinatorial background contribu- 
tions. In the case of the D + — > K* n + v and combinatorial backgrounds, smoothed 
shapes been used. The fit histogram is shown in both plots for reference. 

4 Systematic Studies 

Several studies have been performed in order to assess systematic contribu- 
tions to the uncertainty in the ratio. Three possible contributions have been 
identified. 

The first contribution is due to the final cut selection used to determine the 
branching ratio. In order to estimate a contribution from this source, we vary 
our cuts around the final choice to exercise likely differences between the 
signal and background. Since the D + is longed lived compared to sources of 
background from Df, and other short-lived backgrounds such as those coming 
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from non-charm sources, we vary the significance of the separation between 
the production and the decay vertices from 10a to 20a, and out of target 
requirements for the decay vertex from 0<7 to 2a. To look for poorly formed 
vertices and vertices that are formed from particles that decay into muons 
early in the spectrometer, we vary the confidence level of the secondary vertex 
from 1% to 10%. We have estimated the feed-down from D + — > K~n + fi + u 
using our Monte Carlo simulation, but look for backgrounds we might have 
missed by varying the Cerenkov identification cuts for the pions from 4 to 6 
units of likelihood. The level of the muon misid was checked by changing the 
muon identification confidence level from 1% to 10%, the muon momentum cut 
from 10 GeV/c to 20 GeV/c, and selecting events that left hits in all 6 of the 
muon planes. A very stringent test which dramatically changes the background 
level is to relax the visible mass cut. Though the statistical significance of the 
result suffers due to the inclusion of so much background, this is an important 
check on backgrounds we might have missed coming from higher multiplicity 
modes, which are expected to be small, and combinatorial sources. 

The cut systematic is assessed by measuring the branching ratio for the dif- 
ferent cut combinations and calculating the sample variance for the returned 
values of BR(^A^) _ Because our tested cuts have succeeded in delivering 

BR(D+—*K (t+v) to 

a broad range of signal to background values as well as changes in the final 
yield, this method is likely to deliver a conservative estimate of the systematic 
error due to our cut selection. We find no significant change in the branching 
ratio due to our particular cut choice and assign a systematic uncertainty of 
0.0023 due to our cut selection. 

The second contribution is related to our fit. In order to check the fit against 
possible biases as well as the accuracy of the statistical error reported by 
the fit, we perform the fit multiple times after fluctuating each bin of the data 
histogram using a Poisson distribution. We find that both the mean and width 
of the distribution of the fit results are in agreement with our reported results. 

We test the effect of the fit inputs by changing the D + — > K~7i + fi + u ECY 
used to estimate the amount of background due to K/tt misidentification in 
our D + — > p°fi + u sample by a factor of two and by fitting with no restrictions 
on the D + — > ujfi + u yield. We test the combined Df shape used in the fit 
by fitting our signal using the individual shapes of the modes mentioned 
earlier. In this case, we replace the yield parameter in the fit with a 
parameter representing the D+ — > (f)n + v ECY and we extract the individual 
yields using the branching ratios of these modes relative to Df — >• 0/i + za These 
branching ratios are then varied by ±la. We have also changed the shape of 
the combinatorial background by replacing it with the shape obtained when 
two same sign pions are used to form anirfi vertex. As a final check on the fit, 
we have changed the binning scheme and mass range used in the fit. As in the 
case of cut variations, we calculate the sample variance of the returned values 
and quote this as our systematic contribution. We find this contribution to be 
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0.0038 mostly coming from the uncertainty in the combinatorial background 
shape. 

The third category includes a search for additional, unaccounted for, system- 
atic uncertainty that may come from the detector simulation and/or the charm 
production mechanism. We estimate this by splitting our sample into three 
pairs of statistically independent sub-samples. A powerful test for the produc- 
tion model, trigger, and detector simulation for the signal as well as for the 
background is to split the data according to the D ± momentum. Another test 
of the production model is to look whether we have a decay of a D + or a 
D~ . Our final split sample looks at two different detector configurations. For 
roughly 30% of the FOCUS running we ran without the interleaved silicon 
planes within the target. With this split, we tested the detector simulation as 
well as lifetime dependent backgrounds. 

In order to separate the likely systematic error contribution in any differences 
in the results from larger statistical fluctuations due to reduced statistics, 
FOCUS uses a technique based on the S-factor method of the PDG. In this 
method the branching ratio is measured for each pair of split sub-samples and 
a scaled variance is calculated. The split sample contribution is the difference 
between the scaled variance and the statistical variance if the scaled variance 
is greater than the statistical variance for the entire sample. We find that 
no additional contribution to the systematic uncertainty is indicated by this 
search. 

The total systematic uncertainty is obtained by adding in quadrature all of 
these contributions as summarized in Table 2. 



Systematic Source 


Error 


Cut variations 
Fit variation 
Total 


0.0023 
0.0038 
0.0044 



Table 2 

Sources of systematic errors. The three sources are added in quadrature to obtain 
the total systematic error. 



5 Conclusions 



From 320±44 D+ -> p /2 + u decays and 11, 372±161 D+ -> K'Tr+fi+u decays, 
we report a measurement of the branching ratio 

BR(D + -»• p°u + is) , . . . 
i 4£ — '- = 0.041 ± 0.006 stat. ± 0.004 syst. . 
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Using this result along with the FOCUS measurement of the ratio J7k - J+n +) 
the PDG [16] value of the absolute branching fraction of the decay D + — > 
K~it + it + , and the FOCUS measurement of the D + lifetime [26], we calculate 



T(D + -> p°fi + u) = (0.22 ± 0.03 ± 0.02 ± 0.01) x 10 10 s" 1 



where the last error is a combination of the uncertainties on the quantities 
not measured in this work. When calculating the partial decay width, we have 

corrected ^^+^*k- t!+tt+) w ^ n the updated value for the S-wave non-resonant 
contribution [24]. In Table 3 and Table 4, we compare our result to previous 
experimental results and theoretical predictions, respectively. 

Our result for BR (- p+ ~^ v +v ) [ s consistent with a recent CLEO collaboration 

result on the absolute branching ratios for D + semi-electronic decays [22] and 
represents a significant improvement to the world average for the semi-muonic 
mode. The experimental results indicate that the QCD Sum Rule predictions 
for D + — > p°£ + u [1,5] are too low. 



Reference 


BR(£>+-W+i/) 
BR(B+->FV+i/) 


BR(D+^/e+i/) 
BR(O+-tF e+i') 


E653 [13] 


n 044+ - 034 
u.u^_ 029 




E687 [14] 


0.079 ± 0.023 




E791 [15] 


0.051 ±0.017 


0.045 ±0.017 


CLEO [22] 




0.038 ± 0.008 


This result 


0.041 ± 0.007 





Table 3 

Experimental results for the branching ratio. Statistical and systematic errors have 
been added in quadrature. 
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1-? o to yon oo 




K, 




1 yly > p t u JyL\j h J 


Ball [1] (SR) 


e 




0.06 ± 0.02 


APE [2] (LQCD) 


I 


0.043 ± 0.018 


0.3 ± 0.1 


t r*-»1 / i — V TV r\ 

Jaus [3] (QM) 


i 


0.030 


0.16 


ISGW2 [4] (QM) 


e 


0.023 


0.12 


Yang-Hwang [5] (SR) 


e 


0.018 ± 0.005 


u - u/ -0.02 


O'Donnell-Turan [6] (LF) 




0.025 




Melikhov [7] (QM) 


I 


0.027, 0.024 


0.15, 0.13 


Ligeti-Stewart-Wise [8] 


i 


0.044 




~t 7~ i j i m l ■ To 1 / r t — ^ \ 

Kondratyuk- 1 chem [9] (LI J 


i 


0.035, 0.033, 0.033, 0.032 


0.19, 0.20, 0.18, 0.19 


Melikhov-Stech [10] (QM) 


I 


0.035 


0.21 


Wang- Wu-Znong [1 1J 


o 
i 


0.035 ± 0.011 


0.1 f ± U.U4 


Fajfer-Kamenik [12] 


i 


0.045 


0.25 


This result 




0.041 ± 0.007 


0.22 ± 0.04 



Table 4 



Theoretical predictions for the branching ratio and partial decay width. Most of the 
theoretical predictions are calculated for D° -» p~i + v. To compare these predictions 
with our result, we have used the isospin conjugate relation T(D + — > p°£ + u) = 
1/2 r(D° -» p-^+i/). 
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